Disorder and quantum fluctuations in superconducting films in strong magnetic fields 
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We find that the upper critical field in a two-dimensional disordered superconductor can increase 
essentially at low temperatures. This happens due to the formation of local superconducting islands 
weakly coupled via the Josephson effect. The distribution of the superconducting islands is derived. 
It is shown that the value of the critical field is determined by the interplay of the proximity effect 
and quantum phase fluctuations. We find that the shift of the upper critical field is connected with 
the pinning properties of a superconductor. 

PACS numbers: 74.40. +k, 74.76.-w 



The theory of the upper critical field in homoge- 
neous superconductors was developed long ago 
In the framework of the BCS mean field theory, the 
value of the upper field is weakly temperature dependent 
H c2 (0) - H c2 {T) oc (T/T c0 ) 2 , if T « T c0 . 

However, some experiments (J] show an anomalous up- 
ward curvature of the upper critical field temperature 
dependence as T — > 0. A possible explanation for this 
behavior is based on the following arguments: One can 
not avoid having inhomogeneities in a superconductor. 
The critical field depends on disorder. Local optimal re- 
alizations of disorder may lead to the formation of spatial 
regions where the local upper critical field exceeds the av- 
erage value. These regions form superconducting islands 
weakly coupled via the Josephson effect. At low tem- 
perature, proximity coupling is long-ranged and, thus, 
the global superconductivity may be established in the 
system. For the first time, this idea was proposed by 
Spivak and Zhou H|). They considered mesoscopic ef- 
fects in a dirty superconductor. The islands were formed 
due to the local fluctuations of the density of ordinary 
impurities. It was shown that in the absence of quantum 
fluctuations, even this weak disorder leads to the macro- 
scopic superconductivity at T = for an arbitrarily high 
external magnetic field. 

Similar phenomena at finite temperatures were con- 
sidered in the work of Ioffe and Larkin 0. In this pa- 
per, the mean free path was supposed to be a random 
variable. The optimal fluctuations corresponding to the 
appearance of superconducting islands above the critical 
temperature were found. Thermal fluctuations destroyed 
Josephson coupling. If the energy of the Josephson in- 
teraction between two superconducting islands exceeded 
temperature, than these islands belonged to a same su- 
perconducting cluster. The appearance of an infinite su- 
perconducting cluster corresponded to the genuine tran- 
sition point. 

At low temperatures, the classical thermal fluctua- 
tions are negligible. However, Josephson coupling can be 
suppressed by quantum fluctuations. Fluctuation-driven 
quantum phase transitions in granular superconductors 
were considered in Refs. [0,01. 



In this Letter, we study the suppression of supercon- 
ductivity by the quantum phase fluctuations. We con- 
sider a system of superconducting islands in a strong 
perpendicular magnetic field at low temperatures. We 
examine two different physical mechanisms of disorder 
leading to the formation of the local superconducting is- 
lands. First, we study a dirty two-dimensional super- 
conductor with a random distribution of impurities. We 
show that in the case of weak disorder, mesoscopic fluc- 
tuations may become important only in the very vicinity 
(inside the Ginzburg region) of the old transition point, 
i.e. at H — H C 2(0) ~ H C 2(0)/g, where g — afi/e 2 is the 
dimensionless conductance and the overline means aver- 
aging over disorder hereinafter. We suppose that the sys- 
tem possesses a large dimensionless conductance g >• 1 
and, thus, the Ginzburg region is very narrow. In this 
small region the quantum superconducting fluctuations 
are essential, i.e. the fluctuations of the modulus of the 
order parameter must be taken into account |]i"of . 

Second, we consider a generic disordered system with 
a random diffusion coefficient: 



D(r) = D + 5D(r), 



(1) 



where a short-scale disorder characterized by the Gaus- 
sian white noise is introduced 



SD(r)8D(r') 



D 2 d' 



S(t 



(2) 



The randomness can be connected with localized and ex- 
tended defects present in a superconductor. Deep into 
the superconducting state, this randomness leads to the 
collective pinning effects. Thus, phenomenological con- 
stant d is directly connected with the pinning properties 
of a superconductor JHJ and can be extracted indepen- 
dently from experiments. Let us also note that the model 
can be realized in a system of superconducting grains. Ir- 
regularities of grain sizes and shapes may lead to a ran- 
domness in the tunneling constant between the granules 
and to a random conductance. In this case, constant d 
is of the order of the grain size times a dimensionless 
strength of disorder. In the model defined by Eqs. (Q) 
and (|^), we find the shift of the upper critical field as 
follows: 



1 



(3) 



where Lh = (2eH) is the magnetic length. 

Distribution of superconducting islands. - Let us con- 
sider the vicinity of the BCS upper critical field H C 2(Q) at 
zero temperature. The transition is controlled by the di- 
mensionless parameter h — \H — H C 2(0)] /H C 2(0). The 
fluctuation region in an homogeneous superconductor is 
determined by the condition h < 1/g. We suppose that 
the external magnetic field is such that h < I but lies 
outside the fluctuation region h > g . 

As we have pointed out, even above H C 2 superconduct- 
ing islands appear. To find the distribution of the is- 
lands, we will consider the equation for the order param- 
eter. To obtain this equation, the interaction term in the 
BCS Hamiltonian should be decoupled via a Hubbard- 
Stratonovich field A. Then, the one-electron degrees of 
freedom can be integrated out and one gets an effective 
action for the superconducting order parameter. In the 
vicinity of the transition, an expansion on the order pa- 
rameter is possible and we obtain the following action: 

— Sa = ^ J A*(xi)£~ 1 (xi,X2)A(x 2 )dxidx2 

+ - [ A*(x 1 )A*(x 2 )B({x i })A(x 3 )A(x 4 )]Jdx i , (4) 
J i 

where we use x — (r, t) and dx = d 2 r dt for brevity. 
In Eq.(||), operator £ is the fluctuation propagator: 



r — 



A" 



IT 



A is the BCS interaction con- 



stant, operator IT in the coordinate representation has the 
form II w (r,r') = T^ £ II w (r,r';e) where ir,(r,r';e) = 
C/ e (r, r')(? w _ £ (r, r') and Q e is the Matsubara Green func- 
tion. Let us emphasize that operator II = n+JII consists 
of a mean part and a random part <5LI which is responsible 
for the effects under consideration. Non-linear operator 
B in Eq.(|j) corresponds to the diagrams calculated ex- 
plicitly by Maki |l2| and Caroli et al. |i~3| . In the vicinity 
of the transition we can neglect the randomness in the 
A 4 term. 

The saddle point approximation 5S/5A(r,t) = re- 
sults in the time-dependent Ginzburg-Landau equation 
for a gapless superconductor JlJ|. When considering 
the spatial distribution of the islands we can disre- 
gard dynamic effects and consider the static form of the 
Ginzburg-Landau equation. 

If we neglect randomness in the Cooperon, indeed, 
there are no non-trivial solutions for the corresponding 
mean-field equation above the BCS upper critical field 
h > 0. However, if the random part of kernel 511 pos- 
sesses an eigenvalue which is greater than h a non-trivial 
solution appears. This corresponds to the appearance of 
a local superconducting island. 



To find the distribution of the islands, one should find 
the distribution function of the eigenvalues for the ran- 
dom operator II: 



~ [n{r,r')ip(r')d 2 r' = (e + £e)^(r), 



(5) 



where e = e + 5e is the dimensionless eigenvalue of the 
Cooperon and v is the density of states per spin at the 
Fermi line. In the absence of a randomness in LI, the 
eigenvalues e are well-known. The spectrum is discrete 
and is parameterized by the Landau level indexes. The 
random part smears out the eigenvalues. 
The "density of states" can be defined as 



p(e) = J V {SH} 6(e-e [II]) w [511] 



(6) 



where w [511] is the distribution function for the 
Cooperon which is supposed to be Gaussian with cor- 
relator 5II(ri, r2)<5n*(r3, r4). 

To find the density of states p(e), we use the optimal 



fluctuation method 
functional integral 



1 1 51 . This means that we evaluate 
in the saddle-point approximation. 
Let us note that the problem of finding p(e) is analogous 
to the problem of density of states of a particle in a ran- 
dom potential. In the presence of an external magnetic 
field in 2D, the problem is simplified, since the coordi- 
nate dependence of the wave functions is dictated by the 
magnetic field 0. 

In the case (5e) 2 <C (^e) 2 <§C I, the solution has a form 
of rare islands. In the vicinity of a spherically symmetric 
island located at a point r^, the "wave function" can be 
taken in the following form: 



f 



exp 



(r-r,) 2 



V2ttL h 

In the first approximation we obtain: 



4 7.2 



(7) 



1 If 
5e = -5U QQ = - / ip l (r 1 )5n(r 1 ,r2)ipt(r 2 )d 2 rid 2 r2. 



The distribution function reads: 



p(e) oc exp 



(5ef 
21 



(8) 



where / = <5LIg /z/ 2 . 

In the case of a dirty metal, correlator / can be cal- 
culated with the help of the conventional cross diagram 
technique This yields the following estimate for the 
correlator: I^ 1 ~ g 2 , where g is the dimensionless con- 
ductance and index "1" refers to the first model we con- 
sider (weak mesoscopic fluctuations in a dirty metal). 

In the system with a short-scale randomness in the 
diffusion coefficient (|2|), we can calculate the correlator 
using the differential equation for the Cooperon which in 
the presence of an external magnetic filed has the form: 



2 



[d D(r) d + ie] H(r, r'; e) = (2tw) 5(r - r') 



(9) 



where d = — iV — 2eA(r). One can solve Eq.® using 
a simple perturbation theory with respect to SD. With 
the help of Eq.®, we get the correlator and find the dis- 
tribution function which can be written in the form (^) 
with IjT 1 = 8tt (L H /d) 2 . 

The modulus of the order parameter in a supercon- 
ducting island is random and parameterized by random 
variable e (see EqsjjsJ) and (H)). Using the explicit ex- 
pression (see Ref. pM) for the non-linear operator B in 
Eq.®, one can get the following "mean-field" value of 
the order parameter for a spherically symmetric island: 



D 



|A | = Vi^-^V^- 



(10) 



and the coordinate dependence of the order parameter is 
described by Aj(r) = AoV'i(r), where ipi( T ) is defined in 
Eq. (Q). Let us note that the typical size of a supercon- 
ducting island is Lh- The typical distance between the 
islands is exponentially large R ~ Lh exp [/i 2 /4i] . 

Josephson coupling. - For each realization of disorder, 
there is a fixed spatial distribution of the superconduct- 
ing islands. The interaction Hamiltonian for such a sys- 
tem can be obtained from Eq.® and has the standard 
form 



COS ( 



4>j + Aij), 



(11) 



where Jy is the Josephson energy of the interaction be- 
tween islands i and j and is the phase-shift due to 
the magnetic field. The average value of the Joseph- 
son energy is J(R) oc R~ 2 exp [— R/Lh]- Since the typ- 
ical distance between the islands is exponentially large 
compared to Lh, the average Josephson energy is neg- 
ligible. However, as it was shown in jBj, the variation 
of the Josephson energy decays as a power law only 
J%R) oc R- 4 

The Hamiltonian ( Jil"| ) describes a "frustrated" two- 
dimensional XF-model with random bonds. The frus- 
tration comes both from the Josephson energy which is 
random and from the phase difference due to the mag- 
netic field. At zero temperature such a system should 
show a glassy behavior if there are no effects capable of 
destroying phase coherence between the islands. 

Transition point. - To find the transition point, we 
use action (|j) which describes the dynamics of the su- 
perconducting order parameter. 

Let us present the superconducting order parameter in 
the following form 



A( r> *)=5^ |Aoi|V< (r)e** 



(*) 



(12) 



where ipi is defined in (pi). We consider the islands in 
which the modulus of the order parameter is fixed by the 



static mean-field equations ( |i0| ) and only the phase is al- 
lowed to fluctuate. Finally, we obtain the following action 
describing the system of local superconducting islands: 



S 




dt'rji 



COS [&(*)- &(*')] 1 

(t - t' f E c 



dt 



dt'Jij (t - t') cos [&(t) - </>, [t') + A. K 



The coefficient in the dissipative term is random and con- 
nected with the modulus of the order parameter (Jl^) in 
an island 77, — v | Aoi| 2 /8ireDH. Note that the typical 
value of the coefficient is large i] ~ gh. 

We keep the ui 2 term in the action. As we will see 
below, the effective charging energy appears only as a 
high-frequency cut-off. With the logarithmic accuracy, 
the exact value of E c is not important in our problem. 

Let us integrate out the high-frequency degrees of free- 
dom in the action. First, we consider strong enough mag- 
netic fields so that the network of the superconducting 
islands is very dilute and the average Josephson energy 
is exponentially small. In the domain lu 3> J, only the 
first two terms in action ( [Hf ) are important. In this case, 
the action can be written as a sum of single-island ac- 
tions. The phases in different islands fluctuate indepen- 
dently. With the aid of the single-island action one can 
integrate out the high-frequency phase fluctuations using 
the renormalization group developed by Kosterlitz for a 
spin system with long-range interactions ]l6[ |. Since the 
first term in Eq.(p^) is not Gaussian, coefficient r\ gets 
renormalized when integrating out fast variables. The 
corresponding renormalization group equations are iden- 
tical to the ones derived in pjlql. The solution of the 
RG equation for the renormalized "viscosity" coefficient 
has the following simple form: 



n 1 , E c 

7] w = 77 - — In — 



(14) 



Note that this equation is valid unless rj(u>) becomes of 
the order of unity. This happens at times t c ~ uj^ 1 ~ 
E~ x exp [27r 2 77] . At larger times the phase fluctuates 
rapidly. 

If the Josephson interaction between two islands is 
such that Jijt™ ln < 1, then the Josephson term does 
not affect the phase dynamics at any times and can be 
treated as a small perturbation. However, there is always 
a finite probability of finding a pair of islands for which 
Jyt™ m > 1. In this case, Jy stabilizes the fluctuations of 
the relative phase [(pi(t) — <j)j(t)] and the corresponding 
critical time increases t c ~ E^ 1 exp [2tt 2 [rn + r)j)] . 

Recall, that J 2 oc i?~ 4 and the probability distribu- 
tion for random quantity t c is known and determined by 
Eqs. (§) and ©. Thus, one finds the probability of find- 
ing a pair of strongly correlated superconducting islands: 



(13) 



3 



p 



exp 



-2n 2 g h 



-gi 



(15) 



At large fields, this probability is exponentially small. 
As the external magnetic field decreases, the fraction of 
the strongly correlated superconducting islands increases 
and finite size superconducting clusters are formed. At 
some threshold field, the infinite superconducting cluster 
appears and it corresponds to the macroscopic supercon- 
ductivity. We can estimate the location of the transition 
point as a field at which probability (15) is of the order 
of unity. This yields 



h C 2 = —gi. 



(16) 



Let us mention that result (16) can be ob 



taincd more formally by calculating correlator C — 
J (exp [i(f>j(t) — i(j)j(0)]) s dt. At the transition point, 
the correlator diverges. One can perform the virial ex- 
pansion with respect to the density of islands in C. As in 
the theory of liquids and gases and in the theory of spin- 
glasses with RKKY interactions |l7) , the virial expansion 
can not prove the very existence of the transition. How- 
ever, it determines the transition point if there is one. 
The transition is defined as a point at which all terms 
of the virial expansion become of the same order. Com- 
paring the contribution in correlator C from independent 
islands and the one from pairs, we find Eq. (|l6|). 

In the case of the weak mesoscopic disorder I\ ~ g~ 2 
and the shift of the upper critical field, if any, is small: 
h C 2 ~ 9 • Let us note that this result may acquire 
some logarithmic corrections of the order of (g~^lng) 
which are, however, beyond the scope of our investiga- 
tion. Within the logarithmic accuracy, we can not distin- 
guish the mesoscopic effects under consideration from the 
usual superconducting fluctuations inside the Ginzburg 
fluctuation region. 

In the case of strong disorder, the shift of the crit- 
ical region can be large and it is described by Eq.([|). 
Phenomenological constant d measures the strength of 
disorder which can be connected with dislocation clus- 
ters, grain boundaries in polycrystal samples etc. Let us 
emphasize, that the pinning parameters in a supercon- 
ductor are determined by d. For example, the critical 
current of a superconducting film |l8|| in the collective 



pinning regime is j c /j c o 



H c2 (0) 
H 



d 2 /L 2 Hc2 , where j c0 is 



the depairing current in zero field. Let us also note that 
a possible randomness in the BCS interaction constant A 
would lead to the same effects on the upper critical field. 

At finite temperatures, there is no phase transition in 
a strict sense. At any fields, a finite, though exponen- 
tially small, resistance exists in a two-dimensional super- 
conductor if a magnetic field is applied. One can define 
the upper critical field as a field at which a sharp fall 
in the resistance takes place. At very low temperatures, 



h c i is determined by Eq. (||). As the temperature in- 
creases, h C 2(T) decreases very rapidly. First of all, at a 
finite temperature, the Josephson coupling decays expo- 
nentially at the distances larger than y/ ' D/T. Second, 
the thermal fluctuations destroy the Josephson coupling 
at Jij ~ T . The both effects lead the following estimate 
of the transition temperature T c (h) ~ T c q exp [~h 2 /AI~\ . 
Thus, in a relatively wide region \fl < h < h C 2(0), the 
critical temperature depends on the external magnetic 
field exponentially. 

The increase of H C 2 at low temperatures has been ob- 
served in a number of experiments. The mechanism pro- 
posed in the present paper can give a possible explanation 
for the effect. To reveal whether or not the fluctuation 
effects in disorder are responsible for the increase of H C 2, 
it would be interesting to investigate possible correlations 
between the value of the upper critical field at low tem- 
peratures and the pinning properties. 
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